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ABSTRACT 

We compute the one loop fermion self-energy for massless Dirac + Einstein in 
the presence of a locally de Sitter background. We employ dimensional regu- 
larization and obtain a fully renormalized result by absorbing all divergences 
with BPHZ counterterms. An interesting technical aspect of this compu- 
tation is the need for a noninvariant counterterm owing to the breaking of 
de Sitter invariance by our gauge condition. Our result can be used in the 
quantum-corrected Dirac equation to search for inflation-enhanced quantum 
effects from gravitons, analogous to those which have been found for massless, 
minimally coupled scalars. 
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1 Introduction 

In this paper we compute and renormalize the one loop quantum gravita- 
tional corrections to the self-energy of massless fermions in a locally de Sit- 
ter background. The physical motivation for this exercise is to check for 
graviton analogues of the enhanced quantum effects seen in this background 
for interactions which involve one or more undifferentiated, massless, min- 
imally coupled (MMC) scalars. Those effects are driven by the fact that 
inflation tends to rip virtual, long wavelength scalars out of the vacuum and 
thereby lengthens the time during which they can interact with themselves or 
other particles. Gravitons possess the same crucial property of masslessness 
without classical conformal invariance that is responsible for the inflationary 
production of MMC scalars. One might therefore expect a corresponding 
strengthening of quantum gravitational effects during inflation. 

Of particular interest to us is what happens when a MMC scalar is Yukawa 
coupled to a massless Dirac fermion for non-dynamical gravity. The one loop 
fermion self-energy has been computed for this model and used to solve the 
quantum-corrected Dirac equation 



Powers of the inflationary scale factor a = e^^ play a crucial role in under- 
standing this equation for the Yukawa model and also for what we expect 
from quantum gravity. The Yukawa result for the self-energy consists of 
terms which were originally ultraviolet divergent and which end up, after 
renormalization, carrying the same number of scale factors as the classical 
term. Had the scalar been conformally coupled these would be the only con- 
tributions to the one loop self-energy. However, minimally coupled scalars 
also give contributions due to inflationary particle production. These are 
ultraviolet finite from the beginning and possesses an extra factor of a In (a) 
relative to the classical term. Higher loops can bring more factors of ln(a), 
but no more powers of a, so it is consistent to solve the equation with only 
the one loop corrections. The result is a declining oscillatory behavior in 
the loop-corrected mode functions that seems to betoken the evolution of a 
nonzero fermion mass. A recent one loop computation of the Yukawa scalar 
self-mass-squared indicates that the scalar which catalyzes this process can- 
not develop a large enough mass quickly enough to prevent the super-horizon 
fermion modes from fully experiencing this effect [2]. 
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Analogous graviton effects should be suppressed by the fact that the 
hjj.iyijjip interaction carries a derivative as opposed to the undifferentiated 
Yukawa interaction. What we expect is that the corresponding quantum 
gravitational self-energy will consist of two terms. The most ultraviolet sin- 
gular one will require higher derivative counterterms and will end up, after 
renormalization, possessing one less factor of a than the classical term. The 
less singular term due to inflationary particle production should require only 
lower derivative counterterms and will be enhanced from the classical term 
by a factor of ln(a). This would give a much weaker effect than the anal- 
ogous term in the Yukawa model, but it would still be interesting. And 
note that any such effect from gravitons would he universal, independent of 
assumptions about the existence or couplings of unnaturally light scalars. 

Dirac -|- Einstein is not perturbatively renormalizable [S], however, ul- 
traviolet divergences can always be absorbed in the BPHZ sense j3J IHl El ■ 
A widespread misconception exists that no valid quantum predictions can 
be extracted from such an exercise. This is false: while nonrenormalizability 
does preclude being able to compute everything, that is not the same thing as 
being able to compute nothing. The problem with a nonrenormalizable the- 
ory is that no physical principle fixes the finite parts of the escalating series 
of BPHZ counterterms needed to absorb ultraviolet divergences, order-by- 
order in perturbation theory. Hence any prediction of the theory that can be 
changed by adjusting the finite parts of these counterterms is essentially ar- 
bitrary. However, loops of massless particles make nonlocal contributions to 
the effective action that can never be affected by local counterterms. These 
nonlocal contributions typically dominate the infrared. Further, they cannot 
be affected by whatever modification of ultraviolet physics ultimately results 
in a completely consistent formalism. As long as the eventual fix introduces 
no new massless particles, and does not disturb the low energy couplings 
of the existing ones, the far infrared predictions of a BPHZ-renormalized 
quantum theory will agree with those of its fully consistent descendant. 

It is worthwhile to review the vast body of distinguished work that has 
exploited this fact. The oldest example is the solution of the infrared prob- 
lem in quantum electrodynamics by Bloch and Nordsieck fS| , long before that 
theory's renormalizability was suspected. Weinberg [9 was able to achieve a 
similar resolution for quantum gravity with zero cosmological constant. The 
same principle was at work in the Fermi theory computation of the long range 
force due to loops of massless neutrinos by Feinberg and Sucher (TUIIII]. Mat- 
ter which is not supersymmetric generates nonrenormalizable corrections to 
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the graviton propagator at one loop, but this did not prevent the computa- 
tion of photon, massless neutrino and massless, conformally coupled scalar 
loop corrections to the long range gravitational force More re- 

cently, Donoghue [121 El has touched off a minor industry [THl (HI 1201 HH 122] 
by applying the principles of low energy effective field theory to compute 
graviton corrections to the long range gravitational force. Our analysis ex- 
ploits the power of low energy effective field theory in the same way, differing 
from the previous examples only in the detail that our background geometry 
is locally de Sitter rather than flat.^ 

That summarizes why the exercise we have undertaken is both valid and 
interesting. In the next section we work out the fermionic part of the Feyn- 
man rules for Dirac + Einstein. Section 3 is devoted to the issues associated 
with the graviton propagator. A major complication concerns the impossi- 
bility of employing a de Sitter invariant gauge condition |211 12S] • We give 
a short review of the complex literature on this issue. Then we introduce 
a noninvariant gauge fixing term, isolate the subgroup of de Sitter transfor- 
mations that it respects, and present the gauge-fixed graviton propagator. 
The section closes with a discussion of the BPHZ counterterms necessary for 
our computation. In section 4 we evaluate the contributions from diagrams 
involving a single 4-point interaction. In section 5 we evaluate the more dif- 
ficult contributions which involve two 3-point interactions. Renormalization 
is accomplished in section 6, and our conclusions are given in section 7. 

2 Fermions in Quantum Gravity 

The coupling of gravity to particles with half integer spin is usually accom- 
plished by shifting the fundamental gravitational field variable from the met- 
ric Qfiui^) to the vierbein e^m{.x)? Greek letters stand for coordinate indices 
and Latin letters denote Lorentz indices, and both kinds of indices take val- 
ues in the set {0, 1,2,..., (D — 1)}. One recovers the metric by contracting 
two vierbeins into the Lorentz metric rf"^, 



The coordinate index is raised and lowered with the metric {e^ b = g'^'^^ub), 
while the Lorentz index is raised and lowered with the Lorentz metric (e^ ^ = 

^For another recent example in a nontrivial cosmology see |28| . 
^For another approach see j2(j| . 
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rf"^e^c)- We employ the usual metric-compatible and vierbein compatible 
connections, 

= =^ Af,cd = e''c{e^d,ii - '^^uepd) ■ (4) 

Fermions also require gamma matrices, 7^^. The anti-commutation rela- 
tions, 

{7^7'=}^ (7V + 7V) = -2r^^^J, (5) 

imply that only fully anti-symmetric products of gamma matrices are actually 
independent. The Dirac Lorentz representation matrices are such an anti- 
symmetric product, 

J^^-^(7V-7Y)-^7[¥. (6) 

They can be combined with the spin connection (0)) to form the Dirac co- 
variant derivative operator, 

V, = d,+ '-A,,,r'' . (7) 
Other identities we shall often employ involve ant i- symmetric products, 

= _ ^ V + ^'V - ^'V , (8) 
= ^tNV^ + ^r^'V - V • (9) 

We shall also encounter cases in which one gamma matrix is contracted into 
another through some other combination of gamma matrices, 

7^, = -DI , (10) 

7V76 = p-2)7% (11) 
^y/^^ = Ar]^'^I - {D-4)j^Y , (12) 

7V7V7f> = 2YYl"+iD-A)YYY ■ (13) 
The Lagrangian of massless fermions is, 

/:Dirac = i^e\^\V^i,^ . (14) 
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Because our locally de Sitter background is conformally flat it is useful to 
rescale the vierbein by an arbitrary function of spacetime a(x), 

ep^ = aei,, =^ e^' = a-'e^\ (15) 

Of course this implies a rescaled metric g^y, 

9,u = a'g,. =^ g^'^ = a-'r''- (16) 

The old connections can be expressed as follows in terms of the ones formed 
from the rescaled fields, 

rV = a-'{5^^a, + S^,a,-r''a^g,,)+T''^, (17) 
A^cd = -a~^{e^c'^^d-e^d^t^c)a,u + Af,cd ■ (18) 
We define rescaled fermion fields as follows, 

\1> = a^ip and = a^^p . (19) 

The utility of these definitions stems from the conformal invariance of the 
Dirac Lagrangian, 

/^Dirac = '^e'^.'^'tV^^^ , (20) 

where = d^ + ^A^cdJ"". 

One could follow early computations about flat space background [23 EH] 
in defining the graviton field as a first order perturbation of the (conformally 
rescaled) vierbein. However, so much of gravity involves the vierbein only 
through the metric that it is simpler to instead take the graviton field to be 
a first order perturbation of the conformally rescaled metric, 

dfiu = Vfii^ + i^h^^ with = 16nG . (21) 

We then impose symmetric gauge {e^b = Chp) to fix the local Lorentz gauge 
freedom, and solve for the vierbein in terms of the graviton, 

= (\/^^)/ = + \t^hpb - \^^h^h^b + . . . (22) 

It can be shown that the local Lorentz ghosts decouple in this gauge and one 
can treat the model, at least perturbatively, as if the fundamental variable 
were the metric and the only symmetry were diffeomorphism invariance j2H] ■ 
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At this stage there is no more point in distinguishing between Latin letters for 
local Lorentz indices and Greek letters for vector indices. Other conventions 
are that graviton indices are raised and lowered with the Lorentz metric 
[h^^^ = rj^^hp^, h^^ = rj^fr]'^'^ hpcr) and that the trace of the graviton field is 
h = Tj^^hp^y. We also employ the usual Dirac "slash" notation, 



(23) 



It is straightforward to expand all familiar operators in powers of the 
graviton field, 
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(24) 
(25) 

-g = -Kh + -n~n~ — -K~n'^" tip^ + . . . (26) 
Applying these identities to the conformally rescaled Dirac Lagrangian gives. 



t-pb 

1 + -Kh + -K^h'^ - -K^hP^h^^ + ... 
2 8 4 



Cm..c = ^^^^ + ^ hm^^-hP'^j^id.^-hpp^^^j^r''^ 



8 4 



per 



--hh 



--h" h 



\ih''f,Kp),a + ^h''^hpp^^ 



^^MJP'^^ I + 0{k^) . (27) 



From the first term we see that the rescaled fermion propagator is the same 
as for fiat space. 



iSj 



f-1) fj_ 



(28) 



where the coordinate interval is Ax'^{x; x') 



\x — x 



■'II 2 



- {\r]-r]'\ - i6y. 



We now represent the various interaction terms in (j^Tj) as vertex operators 
acting on the fields. At order k the interactions involve fields, and 
hai3, which we number "1", "2" and "3", respectively. Each of the three 
interactions can be written as some combination Vj"^ of tensors, spinors and 
a derivative operator acting on these fields. For example, the first interaction 
is, 

'^hm^^ = X "^^-^.Kp = X %^,h^i3 . (29) 



# 


Vertex Operator 


# 


Vertex Operator 


1 




5 




2 




6 




3 




7 




4 




8 





Table 1: Vertex operators contracted into '^i^ jhaph^ 



Hence the 3-point vertex operators are, 

V^t = ^t'^h^^ , V^ = -^^^^^ , V^^ = -^{^^^ J^^^) ^d,, . 

(30) 

The order interactions define 4-point vertex operators U'^fj^'^ similarly, for 
example, 

^K'h'm^^ = l^t^r^'^^hij X ^i-^jKphp, = Uttf" X ^^m^Kphp, . (31) 

o o 

The eight 4-point vertex operators are given in Table ^ Note that we do not 
bother to symmetrize upon the identical graviton fields. 



3 Graviton Propagator and Counterterms 

The gravitational Lagrangian of low energy effective field theory is, 

/:Einstem = " (/^-2)A)v^ . (32) 

The symbols G and A stand for Newton's constant and the cosmological 
constant, respectively. The unfamiliar factor oi D — 2 multiplying A makes 
the pure gravity field equations imply R^^j = hg^y in any dimension. The 
symbol R stands for the Ricci scalar where our metric is spacelike and our 
curvature convention is, 

R ^ g^^R,^ = g^'^ {v"^^^^ - V%^^^ + V%„V\^ - . (33) 

Unlike massless fermions, gravity is not conformally invariant. However, it is 
still useful to express it in terms of the rescaled metric and connection 
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-^Einstein = j^i^a""'^ R-2{D - l)a''-^r'' {a,. 'T' 

-(Z}-4)(D-l)a^-^^'^>^a^-(D-2)Aa^|y^ . (34) 



'"P a 



The factors of a which comphcate this expression are the ultimate reason 
there is interesting physics in this model! 

None of the fermionic Feynman rules depended upon the functional form 
of the scale factor a because the Dirac Lagrangian is conformally invariant. 
However, we shall need to fix a in order to work out the graviton propagator 
from the Einstein Lagrangian ()34j). The unique, maximally symmetric solu- 
tion for positive A is known as de Sitter space. In order to regard this as 
a paradigm for inflation we work on a portion of the full de Sitter manifold 
known as the open conformal coordinate patch. The invariant element for 
this is, 

ds"^ = a^i—dif + dx-dx) where a{ri) = ——— , (35) 
^ ^ Ht] 

and the D-dimensional Hubble constant is H = ^A/ (D — l). Note that the 
conformal time t] runs from —oo to zero. For this choice of scale factor we 
can extract a surface term from the invariant Lagrangian and write it in the 
form 



-^Einstein- Surface = {%-i)Ha''~^^r"r''hp„,,Ko + a""-^ ^gr^~g''"r'' 
^ ^\happhfjcj,j ^hf^pphfj^iy -\- ^happhp^yu i hap^^hfjcju . (36) 

Gauge fixing is accomplished as usual by adding a gauge fixing term. 
However, it turns out not to be possible to employ a de Sitter invariant gauge 
for reasons that are not yet completely understood. One can add such a gauge 
fixing term and then use the well-known formalism of Allen and Jacobson 
[5U] to solve for the a fully de Sitter invariant propagator 1211 1121 IM] • 
However, a curious thing happens when one uses the imaginary part of any 
such propagator to infer what ought to be the retarded Green's function of 
classical general relativity on a de Sitter background. The resulting Green's 
function gives a divergent response for a point mass which also fails to obey 
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the linearized invariant Einstein equation j21]! We stress that the various 
propagators really do solve the gauge-fixed, linearized equations with a point 
source. It is the physics which is wrong, not the math. There must be some 
obstacle to adding a de Sitter invariant gauge fixing term in gravity. 

The problem seems to be related to combining constraint equations with 
the causal structure of the de Sitter geometry. Before gauge fixing the con- 
straint equations are elliptic, and they typically generate a nonzero response 
throughout the de Sitter manifold, even in regions which are not future- 
related to the source. Imposing a de Sitter invariant gauge results in hy- 
perbolic equations for which the response is zero in any region that is not 
future-related to the source. This feature of gauge theories on de Sitter space 
was first noted by Penrose in 1963 j^Sl and has since been studied for gravity 
P3] and electromagnetism [3^. 

One consequence of the causality obstacle is that no completely de Sitter 
invariant gauge field propagator can correctly describe even classical physics 
over the entire de Sitter manifold. The confusing point is the extent of the 
region over which the original, gauge invariant field equations are violated. 
For electromagnetism it turns out that a de Sitter invariant gauge can respect 
the gauge invariant equations on the submanifold which is future-directed 
from the source j37j. For gravity there seem to be violations of the Einstein 
equations everywhere [21]. The reason for this difference is not understood. 

Quantum corrections bring new problems when using de Sitter invariant 
gauges. The one loop scalar self-mass-squared has recently been computed 
in two different gauges for scalar quantum electrodynamics [2H]- With each 
gauge the computation was made for charged scalars which are massless, 
minimally coupled and for charged scalars which are massless, conformally 
coupled. What goes wrong is clearest for the conformally coupled scalar, 
which should experience no large de Sitter enhancement over the flat space 
result on account of the conformal flatness of the de Sitter geometry. This 
is indeed the case when one employs the de Sitter breaking gauge that takes 
maximum account of the conformal invariance of electromagnetism in D = 
3 + 1 spacetime dimensions. However, when the computation was done in 
the de Sitter invariant analogue of Feynman gauge the result was on-shell 
singularities! Off shell one-particle-irreducible functions need not agree in 
different gauges but they should agree on shell [10]. In view of its on- 
shell singularities the result in the de Sitter invariant gauge is clearly wrong. 

The nature of the problem may be the apparent inconsistency between de 
Sitter invariance and the manifold's linearization instability. Any propagator 
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gives the response (with a certain boundary condition) to a single point 
source. If the propagator is also de Sitter invariant then this response must be 
valid throughout the full de Sitter manifold. But the linearization instability 
precludes solving the invariant field equations for a single point source on the 
full manifold! This feature of the invariant theory is lost when a de Sitter 
invariant gauge fixing term is simply added to the action so it must be that 
the process of adding it was not legitimate. In striving to attain a propagator 
which is valid everywhere, one invariably obtains a propagator that is not 
valid anywhere! 

Although the pathology has not be identified as well as we should like, 
the procedure for dealing with it does seem to be clear. One can avoid the 
problem either by working on the full manifold with a noncovariant gauge 
condition that preserves the elliptic character of the constraint equations, 
or else by employing a covariant, but not de Sitter invariant gauge on an 
open submanifold [23]. We choose the latter course and employ the following 
analogue of the de Bonder gauge fixing term of flat space, 

Cgf = —a^'-Y'^F^F, , = V'^ih^p,. - \hp.,, + {D-2)Hah^p5^^ . (37) 

Because our gauge condition breaks de Sitter invariance it will be neces- 
sary to contemplate noninvariant counterterms. It is therefore appropriate 
to digress at this point with a description of the various de Sitter symme- 
tries and their effect upon (j37|) . In our D-dimensional conformal coordinate 
system the |D(Z} + 1) de Sitter transformations take the following form: 

1. Spatial translations — comprising {D — 1) transformations. 

V' = V. (38) 
x'' = x' + e\ (39) 

2. Rotations — comprising |(D — !)(/} — 2) transformations. 

V = V, (40) 
x" = R'^x^ . (41) 

3. Dilatation — comprising 1 transformation. 

7]' = kri , (42) 
x'' = kx\ (43) 
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4. Spatial special conformal transformations — comprising (D—l) trans- 
formations. 



77 = 



X 



V 



l-2e-x+\\9Px-x 



X* — 6'^x-x 



i-2e-x+\\ePx-x 



(44) 
(45) 



It is easy to check that our gauge condition respects all of these but the 
spatial special conformal transformations. We will see that the other sym- 
metries impose important restrictions upon the BPHZ counterterms which 
are allowed. 

It is now time to solve for the graviton propagator. Because its space and 
time components are treated differently in our coordinate system and gauge, 
it is useful to have an expression for the purely spatial parts of the Lorentz 
metric and the Kronecker delta, 



V 



ViJ^u + Sl6l and 



Si: - 6X ■ 



(46) 



The quadratic part of £Emstcin + -^^gf can be partially integrated to take the 
form ^h'^^D^J'^hpcr, where the kinetic operator is. 



pa _ 



2(i;-3)Wo'^o 
1 /D-2 



D. 



and the three scalar differential operators are. 

Da = d,[^g^^d, 
Db = 



1 /L)-2 



-9 



D^D- 



(47) 

(48) 
(49) 

(50) 



The graviton propagator in this gauge takes the form of a sum of constant 
index factors times scalar propagators. 



(51) 
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The three scalar propagators invert the various scalar kinetic operators, 

Dj xiAj{x;x') =i6^{x-x) for I = A,B,C, (52) 
and we will presently give explicit expressions for them. The index factors 



are, 



fJ-'^-'- pa 

rpB 

t^v-l- pa 

rnC 

fJ-i^-'- pa 



'^Vfi{pVa)i 



D~3 



Vpi/Vpa ) 



(p'lv)(p'Ja) 1 

2 



(D-2)(D-3) 



(D-3)5X + r7 



I pv 



(53) 
(54) 

(D-3)5X + Vl -(55) 



With these definitions and equation (jK^ for the scalar propagators it is 
straightforward to verify that the graviton propagator (j^Tj) indeed inverts 
the gauge-fixed kinetic operator, 



D 



pa 



pv 



X I 



pa 



A 



5[^5^)i5''{x-x' 



(56) 



The scalar propagators can be expressed in terms of the following function 
of the invariant length £(x; x') between x^ and x'^. 



y{x]x') = 4 sin^ [-Hi{x] x')) = aa'H^Ax\x;x') 



aa H 



{\i] — r]'\ —iSy 



(57) 
(58) 



The most singular term for each case is the propagator for a massless, con- 
formally coupled scalar ji41j . 



iAcfix; x') 



'!)(- 



4\ -^-1 

^ \ 2 ^ 



(59) 



(47r)T ^2 

The A-type propagator obeys the same equation as that of a massless, min- 
imally coupled scalar. It has long been known that no de Sitter invariant 
solution exists |12] • If one elects to break de Sitter invariance while preserv- 
ing homogeneity ()38II39|) and isotropy ()40II41|) — this is known as the "E(3)" 
vacuum — the minimal solution is |44^ 14' 

iA^(x; x') = iAci{x; x') 



Y{D-1) J D r 



2 ' 



(47r)f r(f) \D-AT{D-r 

W^kh r(n+f) 



— TT cot [—Dj+ m[aa j 



2 

T{n + § + l] 



n- 



D 
2 



-2 T{n + 2) 



'-+2\ 



.(60) 
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Note that this solution breaks dilatation invariance (j42H43|) in addition to 
the spatial special conformal invariance (|44ll45p broken by the gauge condi- 
tion. By convoluting naive de Sitter transformations with the compensating 
diffeomorphisms necessary to restore our gauge condition (jHTj) one can show 
that the breaking of dilatation invariance is physical whereas the apparent 
breaking of spatial special conformal invariance is a gauge artifact j^B]- 

The B-type and C-type propagators possess de Sitter invariant (and also 
unique) solutions, 



T{n + 2) V4 

■A , ,s -A . H^-'^^\, ^r(n+D-3)/|/ 



. (SI) 



- 



(4vr)f„t^,r'~''^ r(n + f) V4. 

They can be more compactly, but less usefully, expressed as hypergeometric 
functions [Tflll^ . 



i/\b{x-x') = ^ \ ^ ^ Fi(D-2,l;-;l-f) , (63) 

{Ait} 2 i (yj 2 4 

zAcix;x') = ^ J ^ ) ^f,{D-3,2;-;1-!^) . (64) 

(47r) 2 i (yj Z 

These expressions might seem daunting but they are actually simple to use 
because the infinite sums vanish in D = 4, and each term in these sums 
goes like a positive power of y{x; x'). This means the infinite sums can only 
contribute when multiplied by a divergent term, and even then only a small 
number of terms can contribute. Note also that the i?-type and C-type 
propagators agree with the conformal propagator in D = 4. 

In view of the subtle problems associated with the graviton propagator 
in what seemed to be perfectly valid, de Sitter invariant gauges [24, 25j, it is 
well to review the extensive checks that have been made on the consistency 
of this noninvariant propagator. On the classical level it has been checked 
that the response to a point mass is in perfect agreement with the linearized. 
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de Sitter-Schwarzchild geometry j2Sl- The linearized diffeomorphisms which 
enforce the gauge condition have also been explicitly constructed • Al- 
though a tractable, Z)-dimensional form for the various scalar propagators 
iAi{x]x') was not originally known, some simple identities obeyed by the 
mode functions in their Fourier expansions sufficed to verify the tree order 
Ward identity jlH] • The full, D-dimensional formalism has been used recently 
to compute the graviton 1-point function at one loop order jHOI- The result 
seems to be in qualitative agreement with canonical computations in other 
gauges A D = 3 + l version of the formalism — with regularization 

accomplished by keeping the parameter 5 7^ in the de Sitter length function 
y{x; x') — was used to evaluate the leading late time correction to the 
2- loop 1-point function j^SllSlI- The same technique was used to compute 
the unrenormalized graviton self-energy at one loop order [S3]. An explicit 
check was made that the flat space limit of this quantity agrees with Cap- 
per's result 156^ for the graviton self-energy in the same gauge. The one 
loop Ward identity was also checked in de Sitter background [S3]. Finally, 
the D = 4: formalism was used to compute the two loop contribution from 
a massless, minimally coupled scalar to the 1-graviton function [37]. The 
result was shown to obey an important bound imposed by global conformal 
invariance on the maximum possible late time effect. 

It remains to deal with the local counterterms we must add, order-by- 
order in perturbation theory, to absorb divergences in the sense of BPHZ 
renormalization. The particular counterterms which renormalize the fermion 
self-energy must obviously involve a single if) and a single ip.^ At one loop 
order the superficial degree of divergence of quantum gravitational contribu- 
tions to the fermion self-energy is three, so the necessary counterterms can 
involve zero, one, two or three derivatives. These derivatives can either act 
upon the fermi fields or upon the metric, in which case they must be orga- 
nized into curvatures or derivatives of curvatures. We will first exhaust the 
possible invariant counterterms for a general renormalized fermion mass and 
a general background geometry, and then specialize to the case of zero mass 
in de Sitter background. We close with a discussion of possible noninvariant 
counterterms. 

All one loop corrections from quantum gravity must carry a factor of ~ 

•^Although the Dirac Lagrangian is conformally invariant, the counterterms required to 
renormaUze the fermion self-energy wiU not possess this symmetry because quantum grav- 
ity does not. We must therefore work with the original fields rather than the conformally 
rescaled ones. 



14 



mass~^. There will be additional dimensions associated with derivatives and 
with the various fields, and the balance must be struck using the renormalized 
fermion mass, m. Hence the only invariant counterterm with no derivatives 
has the form, 

n'^m^^tpy/^ . (65) 

With one derivative we can always partially integrate to act upon the ip field, 
so the only invariant counterterm is, 

K'^m^'ipi'P'ipy/—g . 



Two derivatives can either act upon the fermions or else on the metric to 
produce curvatures. We can organize the various possibilities as follows, 

K'^mil){iT/>y''ijj^—g , H?mRil)il)\/—g . (67) 

Three derivatives can be all acted on the fermions, or one on the fermions 



and two in the form of curvatures, or there can be a differentiated curvature, 

1 



K'e^^^R'^- --g^^-R)^^"^tV,^^ , K^e^^R^.tfj^^'tPV^ . (68) 

Because mass is multiplicatively renormalized in dimensional regulariza- 
tion, and because we are dealing with zero mass fermions, counterterms (j65p . 
flMj) and are all unnecessary for our calculation. Although all four coun- 
terterms ()68|) are nonzero and distinct for a general metric background, they 
only affect our fermion self-energy for the special case of de Sitter back- 
ground. For that case i?^,^ = {D — l)H'^g^i,, so the last two counterterms 
vanish. The specialization of the invariant counter-Lagrangian we require to 
de Sitter background is therefore, 

R 



D{D-r 

— > aiK'^^(i^a-^i^a-^ + jj^^—^)i^^> + a2{D-l)DK^H'^^i^^ . (70) 

Here ai and a2 are D-dependent constants which are dimensionless for D = 4. 
The associated vertex operators are, 

Ciij = aiK^(i^a~^i^a^^i^+H'^i^y=aiK'^(a^^i^d'^a~^y, (71) 

C2ij = a2{D-l)DK''HH^ij . (72) 
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Of course Ci is the higher derivative counterterm mentioned in section 1. It 
will renormalize the most singular terms — coming from the iA^ part of the 
graviton propagator — which are unimportant because they are suppressed 
by powers of the scale factor. The other vertex operator, C2, is a sort of 
dimensionful field strength renormalization in de Sitter background. It will 
renormalize the less singular contributions which derive physically from in- 
flationary particle production. 

The one loop fermion self-energy would require no additional countert- 
erms had it been possible to use the background field technique in background 
field gauge j2Hl IHl IHOl EI]- However, the obstacle to using a de Sitter in- 
variant gauge obviously precludes this. We must therefore come to terms 
with the possibility that divergences may arise which require noninvariant 
counterterms. What form can these counterterms take? Applying the BPHZ 
theorem [H El IHl EI to the gauge-fixed theory in de Sitter background implies 
that the relevant counterterms must still consist of k"^ times a spinor differ- 
ential operator with the dimension of mass-cubed, involving no more than 
three derivatives and acting between \1/ and As the only dimensionful 
constant in our problem, powers of H must be used to make up whatever 
dimensions are not supplied by derivatives. 

Because dimensional regularization respects diffeomorphism invariance, it 
is only the gauge fixing term ()37p that permits noninvariant counterterms.'^ 
Conversely, noninvariant counterterms must respect the residual symmetries 
of the gauge condition. Homogeneity (j38ll39p implies that the spinor dif- 
ferential operator cannot depend upon the spatial coordinate x*. Similarly, 
isotropy ()40II41|) requires that any spatial derivative operators di must either 
be contracted into 7* or another spatial derivative. Owing to the identity. 



we can think of all spatial derivatives as contracted into 7*. Although the 

""^One might think that they could come as well from the fact that the vacuum breaks 
de Sitter invariance, but symmetries broken by the vacuum do not introduce new coun- 
terterms |62| . Highly relevant, explicit examples are provided by recent computations for 
a massless, minimally coupled scalar with a quartic self-interaction in the same locally de 
Sitter background used here. The vacuum in this theory also breaks de Sitter invariance 
but noninvariant counterterms fail to arise even at two loop order in either the expec- 
tation value of the stress tensor |44[ I45j or the self-mass-squared It is also relevant 
that the one loop vacuum polarization from (massless, minimally coupled) scalar quantum 
electrodynamics is free of noninvariant counterterms in the same background |64j . 




(73) 
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temporal derivative is not required to be multiplied by 7° we lose nothing by 
doing so provided additional dependence upon 7° is allowed. 

The final residual symmetry is invariance under dilatations fl42H43p . It 
has the crucial consequence that derivative operators can only appear in the 
form a~^df^. In addition the entire counterterm must have an overall factor 
of a, and there can be no other dependence upon rj. So the most general 
counterterm consistent with our gauge condition takes the form, 



ACr 



(74) 



where the spinor function S{b, c) is at most a third order polynomial function 
of its arguments, and it may involve 7° in an arbitrary way. 

Three more principles constrain noninvariant counterterms. The first of 
these principles is that the fermion self-enery involves only odd powers of 
gamma matrices. This follows from the masslessness of our fermion and 
the consequent fact that the fermion propagator and each interaction vertex 
involves only odd numbers of gamma matrices. This principle fixes the de- 
pendence upon 7° and allows us to express the spinor differential operator 
in terms of just ten constants 



+/?2[(a-V'9o)'(a-V5.)] +/53[(a~V5o)(a"V5.)'] + 

+i/7°f/55(a-S°5o)' + f3J{a-''y%){a-W^^)] + (37{a-W^^)A 



Psia-'^'do) + P9{a-Wdi)] + H^Pio}. (75) 



In this expansion, but for the rest of this section only, we define noncommut- 
ing factors within square brackets to be symmetrically ordered, for example, 

+ l{a-V^o){a-W^^){a-y%) + ^(a-W^^){a-'J%f . (76) 

The second principle is that our gauge condition P7jl becomes Poincare 
invariant in the fiat space limit of if 0, where the conformal time is 
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1] = —e ^^/H with t held fixed. In that hmit only the four cubic terms of 
(f7^ survive, 



lim K^H^aS(iHa)-'^^do,iHa)-^Ydi) = n^p^i^^Oof 

H — ^ ^ I 



(77) 



Because the entire theory is Poincare invariant in that limit, these four terms 
must sum to a term proportional to (7'^c)^)^, which implies. 



(78) 



But in that case the four cubic terms sum to give a linear combination of the 
invariant counterterms (j7ip and (j72p . 



+3 (a-i7%)(a-V'9i)' + {a~'fdiY \ = n'^a-'^a-'^ . (79) 



Because we have already counted this combination among the invariant coun- 
terterms it need not be included in S. 

The final simplifying principle is that the fermion self-energy is odd under 
interchange of x'^ and x'^. 



(80) 



This symmetry is trivial at tree order, but not easy to show generally. More- 
over, it isn't a property of individual terms, many of which violate ()8Up . 
However, when everything is summed up the result must obey (j8Up . hence so 
too must the counterterms. This has the immediate consequence of eliminat- 
ing the counterterms with an even number of derivatives: those proportional 
to /35_7 and to Pio- We have already dispensed with /3i_4, which leaves only 
the linear terms, Ps-9- Because one linear combination of these already ap- 
pears in the invariant (j72j) the sole noninvariant counterterm we require is. 



AiZnon = ^Ca^ where C^ij = asK^HH^.j . 



(81) 
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i 




-1 
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-\t^iWAPp]{x]x) ^6^{x-x') 


2 


i/s:2z[-'3A,^](x;x) ^S^{x-x') 


3 




"^Ap^](x; x) Yd" ^^{x-x') 


4 




ex A 


{x; x) 5^{x — x') 


5 






' {x]x')YJ"^ 


6 




a A 


(x; x')YJ'^'' 5^{x-x') 


7 




a A 




8 


{k^ d'f^i[^pAp„]{x] x') 7^7™ 5^{x-x') 



Table 2: Generic 4-point contractions 



4 Contributions from the 4-Point Vertices 

In this section we evaluate the contributions from 4-point vertex operators 
of Tabled The generic diagram topology is depicted in Fig. 1. The analytic 
form is, 



— I 



^4pt 



x; X 



1=1 



af^Apfj 



{x;x') 6^{x-x') 



(82) 




X 



Fig. 1: Contribution from 4-point vertices. 

And the generic contraction for each of the vertex operators in Table ^ is 
given in Table |21 

From an examination of the generic contractions in Table |2l it is apparent 
that we must work out how the three index factors [a/jT^^] which make up the 
graviton propagator contract into i]"^ and r^"''. For the A-type and S-type 
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index factors the various contractions give 

4 



V 



n(3 



rpA 

pa 



V 



n(3 



D-3 

rpB' 



Vpa , V 

= , r/"" 



rpA 

a/l-Lpa 



D 



D-3 



V ' (83) 



For the C-type index factor they are, 

4 



7] 



7] 



ap 



rnC 
Q-P''- pa 



fjiC 
aP-'- pa 



D-3^ „n.n 2 



D-2 



V<x+p_2)p_3) 



Vf3a 



(85) 



On occasion we also require double contractions. For the A-type index factor 
these are, 



rpA 

aP''- pa 



rjlA 

aP-l-pa 



= D{D-l)-2 



D-l 
D-3 



The double contractions of the -B-type and C-type index factors are. 



^»l3^pa 



rpB 

o^pJ-pa 



ap Ba 



^aP^pa 



rpC 

apJ- 



pa 



, r/^^r/ 



rpB 

o^pJ-pa 

rpC 



= 2{D-1), 
{D^-5D + l 



(86) 



(87) 



= 2 



{D-2){D-3) • 



(D-2){D-3) 

TableElwas generated from Table|21by expanding the graviton propagator 
in terms of index factors. 



aP^pa 



(x: X 



rpA 

apJ- 



pa 



iAAix: x') + 



rpB 

aP''' pa 



iAsix: x') + 



apTp„ 



iAc{x; x' 



(89) 

We then perform the relevant contractions using the previous identities. Re- 
lation © was also exploited to simplify the gamma matrix structure. 

From Table El it is apparent that we require the coincidence limits of zero 
or one derivatives acting on each of the scalar propagators. For the A-type 
propagator these are, 

hm tAA[x; x ) = 



lim dfj_iAA{x; x') = 



(4vr)f r(f) 
(4vr)f r(f) 



-7rcot(^D) +21n(a) !>, (90) 



X Ha6°^ . 



(91) 
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('''dT')'^'[ l'^+{''s')^'^Aix;x)5''{x x') 
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B 


n^[i''^^W+i''^')^iAB{x;x) x') 
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lS^^+TD^§)(S^^i^c{x; x) 6^{x-x') 


8 
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«:2^^8(f^?'iA^(x;x')5^(x x') 
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\''\D-2iiD-3)f-iMh'9',] tAcix; x') S^{x-x') 



Table 3: 4-point contribution from each part of the graviton propagator. 
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The analogous coincidence limits for the B-type propagator are actually finite 
in D = 4 dimensions, 

. / ,x H'^-^ T(D - 1) 1 

hm «Ab(x:x') = tt ^^rrr^ x , (92) 

x'-x ^ (47,)f r(f ) D-2 ^ ^ 

lim df,iAB{x;x') = . (93) 

The same is true for the coincidence limits of the C-type propagator, 

. . ,x H^'^T(D-l) 1 , , 

lim^Ac(x;x) = ^z, ^^^^ x m 9Vn ' ^^^^ 
x'^x (47r)~ r(Y) {D-2){D-6) 

lim (9^iAc(x;x') = . (95) 

x'— >x 

Our final result for the 4-point contributions is given in Table lU It was 
obtained from Table El by using the previous coincidence limits. We have also 
always chosen to re-express conformal time derivatives thusly, 

7% = ^-?. (96) 

A final point concerns the fact that the terms in the final column of Table 0] 
do not obey the refiection symmetry. In the next section we will find the 
terms which exactly cancel these. 



5 Contributions from the 3-Point Vertices 

In this section we evaluate the contributions from two 3-point vertex opera- 
tors. The generic diagram topology is depicted in Fig. 2. The analytic form 
is, 



— z 



3pt 



[x: X 



1=1 



al3 
lik ' 



X I 



j=i 



\X\ X 



(97) 




Fig. 2: Contribution from two 3-point vertices. 
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Table 4: Final 4-point contributions. All contributions are multiplied by 
''""VSdI^ . Wedefine A = f cot(^)-ln(a). 
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_i_^2 ^'^JP^^^ls]{x■ x') YJ'"'df.d'u^U\A x') 



Table 5: Generic Contributions from the 3- Point Vertices. 



Because there are three 3-point vertex operators ()30|) . there are nine ver- 
tex products in (j^Tjl . We label each contribution by the numbers on its vertex 
pair, for example, 



I-J 



= iV^\x) X i S {x-x') X iVH^ix') X i 



apApa- 



[X] X 



(98) 



Table El gives the generic reductions, before decomposing the graviton prop- 
agator. Most of these reductions are straightforward but two subtleties de- 
serve mention. First, the Dirac slash of the fermion propagator gives a delta 
function. 



101 



S 



(x; x') = i5 (x — x') . 



This occurs whenever the first vertex is / = 1, for example. 



1-3 



IK 



■r]^^i(^ X i 



S 



IK 



(x; x') X Yr^d' X I 



x; X 



IK 



Yr^5'^{x-x')d'i 



^ A 



x; X 



(99) 

(100) 
(101) 



The second subtlety is that derivatives on external lines must be partially 
integrated back on the entire diagram. This happens whenever the second 
vertex is J=l or J = 2, for example. 



2-2 



2 ' 



S 



IK 



(x;x') X -—Yid'^^t X i 



apApu 



(102) 
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(103) 



In comparing Table Inland Table |21 it will be seen that the 3-point contri- 
butions with 1=1 are closely related to three of the 4-point contributions. 
In fact the [1—1] contribution is —2 times the 4-point contribution with / = 1; 
while [1—2] and [1—3] cancel the 4-point contributions with / = 3 and 1 = 5, 
respectively. Because of this it is convenient to add the 3-point contributions 
with / = 1 to the 4-point contributions from Table EJ 



^4pt _|_ ^3pt 



Km^-^T{D-1) r {D + 1){D-1){D-A) 



(47r)- 

{D~l){D'^-8D^ + 23D-32) 

3{D^-6D + 8) 3 
^4(D-2)2(D-3)2"8 
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D- 



D-3 
2 

D~3 



2{D-3) 



A 



A 



ai/7° 5^(x-x').(104) 



In what follows we will focus on the 3-point contributions with 1 = 2 and 
1 = 3. 



5.1 Conformal Contributions 

The key to achieving a tractable reduction of the diagrams of Fig. 2 is that 
the first term of each of the scalar propagators iAi{x;x') is the conformal 
propagator iAcf(x; x). The sum of the three index factors also gives a simple 
tensor, so it is very efficient to write the graviton propagator in the form, 



[X\ X 



+ E [puTi^]^SAJ{x■,x'), (105) 



I=A,B,C 



where i6Ai{x; x') = iAi{x; x')—iAcf{x; x'). In this subsection we evaluate the 
contribution to using the 3-point vertex operators (IHUj) and the fermion 
propagator PH|) but only the conformal part of the graviton propagator. 



puApfj 



'^Vl,{pVa)u 
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ci\X, X 
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ciyXj Xj 



(106) 
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Table 6: Contractions from the iA^ part of the Graviton Propagator. 



We carry out the reduction in three stages. In the first stage the conformal 
part (jl(J6|) of the graviton propagator is substituted into the generic results 
from TableEland the contractions are performed. We also make use of gamma 
matrix identities such as Q and, 



S 



{x-,x')j^ = {D-2)t[s]{x-,x') and 7.^"^ = -'-{D-1)Y 



(107) 

Finally, we employ relation whenever ^ acts upon the fermion propaga- 
tor. However, we do not at this stage act any other derivatives. The results 
of these reductions are summarized in Table IHl Because the conformal tensor 
factor [apTpl.] contains three distinct terms, and because the factors of 7" J^^^ 
in Table El can contribute different terms with a distinct structure, we have 
sometimes broken up the result for a given vertex pair into parts. These 
parts are distinguished in Table IHl and subsequently by subscripts taken from 
the lower case Latin letters. 
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In the second stage we substitute the fermion and conformal propagators, 



S 



27rf Ax^ 



(108) 



D 



At this stage we take advantage of the curious consequence of the automatic 
subtraction of dimension regularization that any dimension-dependent power 
of zero is discarded, 

hm zAcf(x; x') = and hm d'iAci{x; x') = . (110) 



In the final stage we act the derivatives. These can act upon the conformal 
coordinate separation Ax^ = x^ — x''^, or upon the factor of (aa')^~~ from 
the conformal propagator. We quote separate results for the cases where 
all derivatives act upon the conformal coordinate separation (Table IZj) and 
the case where one or more of the derivatives acts upon the scale factors 
(Table IHl). In the former case the final result must in each case take the form 
of a pure number times the universal factor. 
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'1111 



{aa'y-^-f^'AXf, 

The sum of all terms in Table [7| is, 

-4e^(.;.0 ^ ^r(f )r(f l)(-2Z,^+5Z,-4)(Z,-l)(aar-*^ . 

(112) 

If one simply omits the factor of (aa')^~~ the result is the same as in flat 
space. Although ()112j) is well defined for x'^ ^ x^ we must remember that 
[S](x; x') will be used inside an integral in the quantum-corrected Dirac equa- 
tion (PJ. For that purpose the singularity at x'^ = x^ is cubicly divergent in 
D = 4 dimensions. To renormalize this divergence we extract derivatives 
with respect to the coordinate x'^, which can of course be taken outside the 
integral in (Q) to give a less singular integrand. 



7^Axu \ 1 



Ax2^ 2{D-l) [Ax2^-2 



(113) 



27 



I 


J 


sub 


Coefficient of -^^^ 


2 


1 







2 


2 


a 





2 


2 


b 


-i(D-2)2(D-l) 


2 


3 


a 





2 


3 


b 


\{D-2f{D-l) 


3 


1 




-{D-lf 


3 


2 


a 


UD-l) 


3 


2 


b 


-|(D-2)2(D-1) 


3 


2 


c 


1{D-2){D-1) 


3 


3 


a 





3 


3 


b 


1{2D-3){D-1) 


3 


3 


c 


-\{D-2){D-l) 



Table 7: iAcf terms in wliicli all derivatives act upon Ax^(x; x'). All contri- 
butions are multiplied by ^T{^)T{^-l){aa'Y~^ . 
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(114) 
(115) 



Expression pi5|) is integrable in four dimensions and we could take D = 4 
except for the explicit factor of 1/(D — 4). Of course that is how ultraviolet 
divergences manifest in dimensional regularization. We can segregate the 
divergence on a local term by employing a simple representation for a delta 
function, 
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The final result for Table [7| is 
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{D-2){D-3){D-A) 



'^S'^ix-x') .(118) 



When one or more derivative acts upon the scale factors a bewildering 
variety of spacetime and gamma matrix structures result. For example, the 
[3— 2]b term gives. 



D-2 



s 



IK 



327r^ 

IK^ 



[x;x')d^iAci{x;x' 
I Ax^ 



aa) 2 



{D~2)Ax^ {D-2)Ha5, 



Ax^ 



327r^ V 2 ) 



Ax2^ 



+ 



2Ax^-2 
{D-2)Ha-i^ 



,(119) 



2Ax 



2D-2 



{D -2ya'HAri-f''AXf, {D-1){D- 2)aHAri-f''Ax^ 



2Ax 



2D 



Ax2^ 
(D-2)2aa'if2^^Aa;^ 

4Ax2«-2 



..(120) 
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rj« 

_n7 


ri LST} 'y^ iXXij. 








2 


1 













2 


2 


a 











2 


2 


b 





-|(i^-2)a' 


UD-2)Da' 


2 


3 


a 











2 


3 


b 





\{D-2)a' 


-l{D-2)Da' 


3 


1 




i(^-l) 








3 


2 


a 


_ 1 








3 


2 


b 


-1(^-2)^ 


i(^-2)« 


i(D-2)V 
-i(D-2)(D-l)a 


3 


2 


c 


-1(^-2) 








3 


3 


a 


l^P-2f 





l(D-2)^(a-a') 


3 


3 


b 


-|(2D-3) 








3 


3 


c 


1^(^-2) 









Table 8: iA^ terms in which some derivatives act upon scale factors. All 
contributions are multiplied by ^^ro r^('f )(aa')^~^. 
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The first term of (|12Uj) originates from both derivatives acting on the con- 
formal coordinate separation. It belongs in Table [7| The next three terms 
come from a single derivative acting on a scale factor, and the final term in 
1)1201) derives from both derivatives acting upon scale factors. These last four 
terms belong in Table |H1 They can be expressed as dimensionless functions 
oi D, a and a' times three basic terms, 



IK 



IGtt 



D 



aa'Y 2 



-{D-2f X '''' 



+ 



l(D-2)V-i(D-l)(Z}-2)a 



X 



121) 



These three terms turn out to be all we need, although intermediate 
expressions sometimes show other kinds. An example is the [3—1] term. 



1 /D-1 



K^d' \ ^iAci{x;x)i S {x\x')'^^ 



8^ 



D-2 



d' 



Ax^ ^2Ax^-2 



7^/3, 

Ax^ 



-r 



{D-2) Ax2^ 2^ ' Ax^D-^ 



Aa;2^-2 



122) 



>.(123) 



As before, the first term in ()123|) belongs in Table [7| The second and third 
terms are of a type we encountered in p20|) but the final term is not. How- 
ever, it is simple to bring this term to standard form by anti-commuting the 
through either 7°, 



aa'H^j^Y^x^j^ = -aa'H^-f^'Ax^ - 2aa'H^Ar]j^ , (124) 
= -aa'H^Y^Xf,-2{a-a)H-f^ . (125) 



Note our use of the identity (a — a') = aa'HArj. 

When all terr 
can be extracted 



When all terms in Table |H1 are summed it emerges that a factor of H^aa' 



IK 



IGtt 



D 



— 



(aa'V-^- 



1 „^ „x aa 

^ 7D + 8)x — 
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--(D-2)(a-a') x^-^--{D-2){3D-2)ia-a') x 



Note the fact that this expression is odd under interchange of and x'^. 
Although individual contributions to the last two columns of Tableware not 
odd under interchange, their sum always produces a factor of a—a' = aa'HAr) 
which makes (jl27|) odd. 

Expression ()127|) can be simplified using the differential identities, 

Ar]'^YAXf, _ dl / Y^x^ 



Ax^D 4(D-2)(D-1) VAx2^-4 



1 'j'^Ax^, ^ 1 7°Ar/ 



2(D-1) Ax2^-2 D~l Aa;2^-2 ' 
7°Ar/ 7°9o / 1 



Ax2^-2 2(D-2) VAx2^- 
The result is, 



+ A^9n-9. ' (128) 

(129) 



—I 



m = 'jfl^rH-) (aa'^'-^ [ p^-llD^ + 23D-12) 7^ 



167r^ V2/' ' [ 16(D-1) Ax2^-2 



D 1 \ 1 /3D-2^^2/7^Aa;, 



16(1^-1) ' "VAx2^-4 

We now exploit partial integration identities of the same type as those pre- 
viously used for Tabled 

^'^^^ = ( ^ \ = J-(J-) + 0(D-4) (131) 

Ax2^-4 2(D-3) VAa;2^-6^ 2 V A,t2^ ^ "^^"^ ' ^'""'^ 

7^Ax^ ^ -^d^ / 1 X 

Ax2^-2 4(D-2)(D-3)(D-4) VAx2^-6J ' 

_^a2rin(^i2Ax^l mf/.^-^ ^^^(x-xQ 

- "1^1 Ax2 / + ^^^"^^ - ^TXfT (/^-3)(/^-4) ' ^^^^^ 

1 92/1 



Ax2^-4 2(D-3)(D-4) VAx2^-6 
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" 4\ Ax2 1+'^^'^ r(f-l) (D-3)p-4) • 

It is also useful to convert temporal derivatives to spatial ones using, 

7°(9o =^ and 
Substituting these relations in ()130|1 gives, 



(133) 



(134) 



^ i -fD3-13D^ + 27L'-12)^ 

2Ht{D-1){D-3){D-A)\ ^ ^ 



-2D(D-2)^U^(x-x') + 



29-3-7r4' 



/ln(/i^Ax' 



+5^(V2-a2)(^)|+0(D-4) . (135) 



5.2 Sub-Leading Contributions from i5A 



A 



In this subsection we work out the contribution from substituting the residual 
A-type part of the graviton propagator in Table El 

2 



a/3 per 



[x: X 



VapVcrl3~^VacrVpf3' 



pa 



Z(5Aa(x; X 



(136) 



T/3'' -laa-ipp £)_3' 

As with the conformal contributions of the previous section we first make the 
requisite contractions and then act the derivatives. The result of this first step 
is summarized in Table IHl We have sometimes broken the result for a single 
vertex pair into as many as five terms because the three different tensors 
in ()136|) can make distinct contributions, and because distinct contributions 
also come from breaking up factors of 7"J^'^. These distinct contributions 
are labeled by subscripts a, b, c, etc. We have tried to arrange them so 
that terms closer to the beginning of the alphabet have fewer purely spatial 
derivatives. 

The next step is to act the derivatives and it is of course necessary to 
have an expression for i6AA{x; x') at this stage. From ()60p one can infer, 

iSA^ix; x') = 

T(^ + l)(aa'?-f H^'^r(D-l)i /vr x ,^ 

-TT cot ( —D 1 + in(aa ) 



_r(f+i) 



Ax 



D-4 



(47r)f 



H^-'^ ^jlT{n+D-l) 



(47r)f^Jn r(n+2 



(4 



n 



r(f) [ ■■--^V2- 

1 r(n+f + l)^yy 



" 2 " 



D 
2 



-2 r n + 2 



(137) 
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sub 


iVj ^{x) i[S\{x; x') iV^ {x') [ai3T^\ ioi^aKx] x') 


2 


1 




-T^)^ 9^§i[S] {x- x')-f^'l5^A{x■, x')} 


2 


2 


a 




2 


2 


b 


+ \K^di{-ikOl>l[S\[X] x')-ikiol\A[x] X )} 


2 


2 


c 




2 


3 


a 


2(D-3) ^^^^['5'](x; x')(JI'i5AA{,x;x') 


2 


3 


b 


-\K^lkdti[S]{x] x')'y(kdi)iSAA{x; x') 


2 


3 


c 


+4(D-3) ^^H'S^J )^«(iAA(x;a;') 


3 


1 


a 


UM)'^ '-W'^^^^^i-'-- )'[5](.r: .rOV'} 


3 


1 


b 


- 2(IJ33) '^'^i.f x')i[S] {x; x')r} 


3 


2 


a 


2(D-3) K'^dki^iSAA{x; x')i[S] (x; x')7fc} 


3 


2 


b 


4(D^_3) K^<9fc{^?5AA(a;; x')i[S]{x; x')7fc} 


3 


2 


c 


^{i [5] (x; x') ^i5A^ (s; x') } 


3 


2 


d 


+|«^<9fc{7^i[S'](a;; x')7A^5Aa(x; x')} 


3 


3 


a 


-\i^3)^'r^[S]{x; x')d,^'tSAA{x; x') 


3 


3 


b 


~4(D-3) i^'^l^^S] [x] x')d^^i6AA{x; x') 


3 


3 


c 


+ 4(D-3) '^'^^k^[S] {x; x')dk^ i5Aa{x; x') 


3 


3 


(1 


-UM)'''^^iAS](r:xy)iJi/)AA{.nx') 


3 


3 


e 


-^K^jki[S]{x; x')jk^HSAA{x; x') 



Table 9: Contractions from the iSAA part of the graviton propagator 
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Function 


Vertex Pair 2-1 


Vertex Pair 2-2 




(D-l) 
(D-2)(D~3)2(D~4) 









-1 




(D-2)(D-3)2(D-4) 


(D-2)(D-3)2(D-4) 




-2 





A,V'^{^;^) 





D(D2-3D-2) 
4{D-2)(D-3)2(D-4) 







(D2-3D-2) 
2(D-3) 







(D-2)(D-3)2 







-2(§5|) 



Table 10: z5A^ terms giving both powers of Ax^. The two coefficients are 
A, ^ ^r(f +i)r(f )(aa')2-^ and ^ $^r(D-2)[in(aa')-.cot(f^)]. 



In D = A the most singular contributions to ()97|) have the form, i^A^i/Ax^. 
Because the infinite series terms in fjl37|) go like positive powers of Ax^ 
these terms make integrable contributions to the quantum-corrected Dirac 
equation (Q). We can therefore take D = A for those terms, at which point 
all the infinite series terms drop. Hence it is only necessary to keep the first 
line of p37|) and that is all we shall ever use. 

The contributions from i5^A are more complicated than those from zA^f 
for several reasons. The fact that there is a second series in (jl37|) occasions 
our Table ITUl These contributions are distinguished by all derivatives acting 
upon the conformal coordinate separation and by both series making nonzero 
contributions. Because these terms are special we shall explicitly carry out 
the reduction of the 2—2 contribution. All three 2—2 contractions on Table IHl 
can be expressed as a certain tensor contracted into a generic form, 

2 -1 r 

5ij5ki + 5ik5jt-^^5ii5jk\ X —di-ij\i5/\A{x]x')dki[S]{x]x')-ii^ . (138) 

So we may as well work out the generic term and then do the contractions 
at the end. Substituting the fermion propagator brings this generic term to 
the form. 

Generic = —di^jli5AA{x; x')dki[S]{x; x')^e\ , (139) 
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= -^^^^5.7,{^5A^(^; ^05.(^^)7.}. (140) 

Now recall that there are two sorts of terms in the only part of i(5(A^(x; x') 
that can make a nonzero contribution for D — A, 

TA / A Ti§ + l){aa'f-T 

z5Aa,(x;x') ^ ^g^^ Ax^-^ ' ^^^^^ 

i5A,,,{x;x') ^ ^^^^{-ncot{^D)+ln{aa')}. (142) 
(47r) 2 i 1,-2 J / 

Because all the derivatives are spatial we can pass the scale factors outside 
to obtain, 



Generic 1 

,^2^2r(|)r(|+i) n r^i . /T^A X 

Generic2 



2^+2 



-^r(L>-l){-7rcot(|i?)+ln(aa')}9,7A(^^)7£ , (144) 



^|^r(D-2){-7rcot(|D)+ln(aa')}a.7,a.^7.(^) ■ (145) 
To complete the reduction of the first generic term we note, 

1 (YAx^\ ^ Ik DY^x^Axk . . 



2VD-2>'Ax20-4^2(L>-2) HAx2^-4/ ' ^-"^^^ 

Hence the first generic term is. 

Generic, ^l^^-lM±iL(„,a')-f 

2%^ (L)-4)(L)-3)(D-2)^ ^ 

X {L'a,7,afe ^7^ - d%jjjkli} . (149) 
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Now we contract the tensor prefactor of (|138|) into the appropriate spinor- 
differential operators. For the first generic term this is, 



D-3 



2 



+ 



■d^^ai^. (150) 



^-3 

This term can be simphfied using the identities, 

f^f = 2^V2 = VV- 2^V^ = -VV+2VV5o , (151) 

= -7.7.^-2^= (1^-1)^-2^= (D-3)^ + 27%, (152) 

= -{D-l)^=_l^l^^ . _ (153) 

lifli = -l^l^p-2f={D-3)f. (154) 

Applying these identities gives, 
2 



D-3 



X 



{Ddi-fjdk ^7£ - d'^di-fj-fkle} 



2D 



D-A 
D-3 



D-3 

For the second generic term the relevant contraction is, 

2 



D-3 



(9V -(155) 



SijSke + SikSje- 



D-3 



D-5^^^^ 



^D-3^ 
D- 



(156) 



Z.-3AV-4(§^)V^^. (157) 

In summing the contributions from Table ^1 it is best to take advantage 
of cancellations between Ai and A2 terms. These occur between the 2nd and 
3rd terms in the second column, the 4th and 5th terms of the 3rd column, 
and the 6th and 7th terms of the 3rd column. In each of these cases the 
result is finite; and it actually vanishes in the final case! Only the first term 
of column 2 and the 2nd term of column 3 contribute divergences. The result 
for the three contributions from [2—1] in Table ITUl is. 

K^H^fi^-' (D-i)r(f +1) 



25vrf {D-3y{D-A) 



{aa'Y--^6''{x-x') 



+ 



2H^ ' 



2 ^ 



|-ln(yU^Aa;^ 
Ax2 



A + 2ln{lH'^Ax^ 



-0(D-4).(158) 
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3 


c 





(Si) 
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2D 


3 
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4(D-l)(D-2) 











3 


1 


b 


2(§z|) 










3 


2 


a 





4(§ei) 








3 


2 


b 




(Si) 


2(Si) 


-4(Si) 


3 


2 


c 




-i(^+l) 




2{d-i) 


3 


2 


d 








2(d-i) 


3 


3 


a 


n(D-l){D-2) 
^ {D-3) 











3 


3 


1) 




rD-U 








3 


3 


c 


VD-3/ 


VD-3/ 








3 


3 


d 


(r,-i)(r.-5) 




(D-5)(D-2) 


(D-5)(D-2) 


4(D-3) 




4(D-3) 


2(Z)-3) 


3 


3 


e 


-1(^-1)^ 


1(^-0 


i(«-2)(i^-i) 


-1(^-2) 



Table 11: iSA^ terms in which all derivatives act upon Ax'^{x;x'). All con- 
tributions are multiplied by ^^r(^+i)r(^)(aa')^""?. 



38 



The result for the five contributions from [2—2] in Table ^1 is, 

K'^H^fi^-^ r(f+i) 



+ 



2H^ 





pln(/i2Aa;2)^ 











2 + ln(i/72^x2) 



Ax2 



0(D-4). (159) 



As might be expected from the similarities in their reductions, these two 
terms combine together nicely in the total for Table EH 



-! lintel' 



/€2i/2^«-4r(f + l)(aa')2-^ 



2%T 



(D-3)2(D-4) 



ln(/i2Ax2) 



Ax2 



/ ^ 9\ r2+ln(ii/2Aa;2^ 



Aa;2 



0(^-4). (160) 



The next class is comprised of terms in which only the first series of iSAa 
makes a nonzero contribution when all derivatives act upon the conformal 
coordinate separation. The results for this class of terms are summarized 
in Table ^2 In reducing these terms the following derivatives occur many 
times, 



diiS^A^x; x) 
9oi5A^(x; x ) 



H 



Ax' 



D 



+ 1 ) (aa 



i\2- 



9qz5A^(x; x) 



STrf 



Ax^-2 
\Ax^ 



-d[i5AA{x;x') (161) 
aH 1 



-2 2Ax^-4 J 

iy^-2r(D-i) 



2^7rT r(f 



Ar; 



-) 



aif, (162) 



Aa;^-2 2Ax^-4 
H^'^V{D-l) 



+ 



2^7rf r(f) 

We also make use of a number of gamma matrix identities, 
7'^7m = -D and = -(^-1) , 



a'i/ . (163) 



(164) 
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{D-2)Y and YYY = {D-1)Y - 2^" , (165) 
-Ax^ and (YAx*)^ = -||Aff , (166) 
(D-l)7^'Ax^ - 27* Ax* , (167) 
||Af|| VAx/, - 2\\AxfYAx' . (168) 



In summing the many terms of Table ITD the constant K = D — 2/{D—3) 
occurs suspiciously often, 



x: X 



^r(y+i)r(-)(aa')-- 



X • 



[-2{D-iy 

'D-2 



D-1 



K 



l^Ax^ 



+ 



-(D-2)- 



K 



YAxi 

Ax2^-2 



|Af|| V^a;^ (D-2)(D-4) \\AxfYAxi 



Ax2^ 



(D-3) 



The last two terms can be reduced using the identities, 



|Ax||27^Ax„ 1 7''Ax 



1 7* Ax, 



Ax2^ 



/ 7'^Ax, 



■.(169) 



Ax2^ 2 Ax2^-2 ^ D-1 Ax2^-2 ^ 4(D - 2) (D - 1 ) V Ax^^-^ 



Ax||27*Axi I/L' + In 7' Ax 



Ax2« 



2^D-l) 



+ 



7*Axi 

Ax2^-2 ' 4(D-2)(Z}-1) VAx2^-4 



,(170) 
(171) 



Substituting these in (jl69j) gives. 



x; X 



2H^ 



■ + l)r(^)(aa') 



2-^ 



y^Ax^ 

Ax2^-2 ^ 



(D-2)(D2_5/} + io) 



2(D-1)4 
DK 1 7'Ax 



2(D-l)(L)-3) 4(D-1) 
KV^ / 7''Ax^ ^ , {D-A)V^ 



16(D-1) VAx 



r.2D-A 



+ 



Ax2^-2 

7''Axi 



4(D-l)(D-3) Ax2^- 



.. (172) 



We then apply the same formalism as in the previous sub-section to partially 
integrate, extract the local divergences and take D = 4 for the remaining, 
integrable and ultraviolet finite nonlocal terms. 
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(D-l)(D-4) 
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3 


3 


c 











Table 12: i5A^ terms in which some derivatives act upon the scale factors of 
the first series. The factor ^-^^V{^+i)V[^){aa'Y~^ multiplies all contribu- 
tions. 



277r# 



P-3)p-4) 
DfsT (L)-2)(D2-5L> + 10) 



4(D-1) 2(L)-l)(D-3) 



29-3-7r4 



-15^9^-4^9^ 



ln(/i^Ax^ 
Ax2 



The final class is comprised of terms in which one or more derivatives act 
upon a scale factor. Within this class we report contributions from the first 
series in Table IT^ and contributions from the second series in Table ITSl Each 
nonzero entry in the 4th and 5th columns of Table IT^ diverges logarithmically 
like l/Aa;^^"'^. However, the sum in each case results in an additional factor 
of a — a' = aa'HA'r] which makes the contribution from Table IT^ integrable. 



€2j 



2HD 



D ^ fD 

2+1)^(2 



aa 



/\3- 



-In 7°A77 
1,) Ax2^-4 



2 V D-?, 



1 /3L'-4n yAxi7^Ax^7°Ar7 {D-l){D-A) 7^Ax^ 



Ax2^-2 



4(D-3) Ax2^-4 



..(174) 
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a 








(D-l)a 
(D-2)(D-3) 


3 


2 


a 










Table 13: idA^ terms in which some derivatives act upon the scale factors of 
the second series. All contributions are multiplied by '^^0+2^0 ^{d-i). 



This is another example of the fact that the self-energy is odd under inter- 
change of and x'^. 

The same thing happens with the contribution from Table 



+ V:^J +^ VD^Jr(f) Ha j-^^^^ 

We can therefore set D = A, at which point the two Tables cancel except for 
the delta function term, 

(176) 

It is worth commenting that this term violates the reflection symmetry ()80|) . 
In D = A it cancels the similar term in ()1()4|1 . 

5.3 Sub-Leading Contributions from iSAb 

In this subsection we work out the contribution from substituting the residual 
-B-type part of the graviton propagator in Table 



af^Apcr 



+ OTlp. + 5l^5"^«p + Wp'lo.. ^S^B . (177) 



As in the two previous sub-sections we first make the requisite contractions 
and then act the derivatives. The result of this first step is summarized in 
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Tfc^oH'S^Jia^; a;')7'' dQidAB{x; x') 


2 


3 


b 




2 


3 


c 


1 2 fco ' TA / /NO • r Ol / /\ 

-^K'^YdktoABix; X ) dQt[S\{x; x )x 
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7 a zpj(x; x j7 Okto/\B[x, x ) 
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o 
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gfi;^c'o|7'''46J(x; x')7fc dozdAij^x; x )| 


o 
o 


2 


b 


gK 7 (^^{^[dJ^x; x j7" aoidAfi^x; X j| 


3 


2 


c 


-i«V5o{^['5](x; x')7' 4^5Ab(x; x')} 


3 


2 


d 


-|'t^5fc{7°^['^](a;; a:')7° <9^«5Ab(x; x')} 


3 


3 


a 


- j^fi;2^fc«[5'] (x; x')-f''dodQi6AB{x; x') 


3 


3 


b 


^K'^'y^i[S]{x; x')7'' Okd^idABix; x') 


3 


3 


c 


-j^K^-i''i[S]{x- x')7° 5o<9fc^5AB(x; x') 


3 


3 


d 


^fi:V^['5](x;x')7°V2z5AB(x;x') 



Table 14: Contractions from the id^B part of the graviton propagator. 



Table EI We have sometimes broken the result for a single vertex pair into 
parts because the four different tensors in ()177j) can make distinct contribu- 
tions, and because distinct contributions also come from breaking up factors 
of j"J^^. These distinct contributions are labeled by subscripts a, b, c, etc. 

the residual of the -B-type propagator ()61|) after the confor- 
mal contribution has been subtracted, 

^5AB[x^x)- ^^^^ ^^^^^ ^^^^ 

H^-^ - rr(n+|).ysn-f+2 T{n + D-2)ryy\ 
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T 
1 


J 


sub 






||Ax||27MAa::,j 


||Ax|p7'Aa;i 






Aa;2D 


Aa;2-D 


2 


3 


a 




-(D+1) 


-D(D-l) 


2D 


2 


3 


b 


(D-1) 


-2D+1 


-D 


2D 


2 


3 


c 





-(D-1) 


-D 


2D 


2 


3 


d 





-1 


-D 


2D 


3 


2 


a 


-2(D-l)(D-2) 


3D-5 


2(D-1)2 


-4(D-1) 


o 
O 


o 
Z 


b 


-(D-1) 


3(D-1) 


2(D-1) 


-4(D-1) 


Q 
O 


9 


c 





2D— 3 


2(D-1) 


-4(D-1) 


3 


2 


d 


-(D-1) 


2D-1 


2(D-1) 


-4(D-1) 


3 


3 


a 


i(D-l)(D-3) 


-(D-3) 


-i(D-l)(D-2) 


(D-2) 


3 


3 


b 







-1(15-2) 


(D-2) 


3 


3 


c 





-i(D-2) 


-1(15-2) 


(D-2) 


3 


3 


d 


i(D-l) 


-(D-1) 


-i(D-2) 


(D-2) 


Total 






-i(D-l)(D-2) 


3(D-2) 


i(D+2)(D-2) 


-4(D-2) 



Table 15: iblS^B terms in which all derivatives act upon Ax^(x;a;'). All con- 
tributions are multiplied by ^^r2(y)(L)-4)(aa')2~f . 



As was the case for the z5AA(a;; x') contributions considered in the previous 
sub-section, this diagram is not sufficiently singular for the infinite series 
terms from i5A5(x;x') to make a nonzero contribution in the D = 4 limit. 
Unlike «5A^(a;; x'\ even the ?t, = terms of z5As(x; x') vanish for D = A. This 
means they can only contribute when multiplied by a divergence. 

Contributions from the [2—2] vertex pair require special treatment to take 
advantage of the cancellation between the two series. We will work out the 
"a" term from Table IT^ 



2-2 




A key identity for reducing the [2 — 2] terms involves commuting two deriva- 
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tives through 1/Ax 



7d„d, 



D-4 



Ax^^- 



(181) 



This can be used to extract the derivatives from the first term of iSAB{x] x'), 
at which point the result is integrable and we can take = 4, 



2-2 



1 iK^H'^D\„^D 



2%^ 



[aa 



^ Ax^-"^ 



(-39o^+7°V2+(i^-l)^fi'o) 



AxD- 



D-3 

X 



-{aa'Y 



do -I {D-A)Ha') (-36*0 ^+7° + {D-l) ^do) 



' 2%^ 



-7 



Ax^ 



+ 0(D-4) . 



(182) 

(183) 
(184) 



Of course the second term of iSAs is constant so the derivatives are already 
extracted, 



2-2 



Ax^- 



2%4 



■7%(3ao' + V' 



Ax^ 



+ 0{D-A) . 



,(185) 
(186) 



Hence the total for [2— 2]^ is zero in D = 4 dimensions! 

The analogous result for the initial reduction of the other [2—2] term is. 



2-2 



m^r(f-l) 



xaJ^(5AB(x;xO(-7"5o5fc^4+7feao'+7fc^7%) 



Ax^- 



(187) 



The results for each of the two terms of iSAb are. 
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2-2 



2%^ 



D-3 



{aa'y 



X 



2-2 



ik'H''-^T{D-2 



(-270aoV2+^V'+^(9i 
1 



1 



2%- 



2^+%^ D-2 



Ax2 

On V72 av72 



27"9oV"-^V^ 



+ 0(1^-4) 
1 



2^4 



1 



Ax2 



0(D-4) . 



, (188) 
(189) 
(190) 
(191) 



Hence the entire contribution from [2—2] vanishes in Z) = 4. 

The lower vertex pairs all involve at least one derivative of iSAs, 



dii6AB{x; x ) 
doiSAsix; x ) 
d'oiSABix] x') 



^'^P(D-4)(aa')-T 



Ax* 



(D-4)(aa' 



167r~ 
lovr 2 



Ax^-2 

\ Ax^-2 2 Ax^-4 
f Ar] a'H 



-d[5AB{x-x') ,(192) 



Ax^-2 2Ax^-4 



(193) 
(194) 



These reductions are very similar to those of the analogous iSAa terms. We 
make use of the same gamma matrix identities p64lll68|l that were used in 
the previous sub-section. The only really new feature is that one sometimes 
encounters factors of Arf which we always resolve as, 



Ar^^ 



-Ax^ + IIAfl 



(195) 



Table IT31 gives our results for the most singular contributions, those in which 
all derivatives act upon the conformal coordinate separation Ax^. 

The only really unexpected thing about Table is the overall factor of 
{D — 2) common to each of the four sums. 



-![E'IISj(rc;i') 



.^r^(f)(Z>-2)(Z>-4)(aaO= 

7*Axj 1 ||Axp7^Ax^ 
'ax2^-2 ^2*^ ^ ' A^ 



\ 2^ ' Ax^D-^ 
II Ax|p7*Ax,- 1 , 
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As with the resuh of Table [TTl we use the differential identities (|17(JH17l() to 
prepare the last two terms for partial integration, 



2 



8(D-l)(D-2) VAx2^-4/ (D-l)(D-2) VAx2^-4^ 
^^(2)^"") ili^^)^^ -8 (D-l)(D-3) 



(197) 



1 (D + 2)(D-4) 1 (Z^-4) ] 1 

"l6(D-l)(D-3)^^ +2(D-l)(D-3)^^ JA^- ^^^^^ 

The expression is now integrable so we can take D = A, 

= !|!|!{_l^#}^ + 0(fl-4) . (199) 

Unlike the iSA^ terms there is no net contribution when one or more 
of the derivatives acts upon a scale factor. If both derivatives act on scale 
factors the result is integrable in D = 4 dimensions, and vanishes owing to 
the factor of (Z) — 4)^ from differentiating both a'^~~ and a''^~~ . If a single 
derivative acts upon a scale factor, the result is a factor of either (D — 4)a 
or [D—A)a' times a term which is logarithmically divergent and even under 
interchange of and x'^. As we have by now seen many times, the sum 
of all such terms contrives to obey reflection symmetry (j8(Jj) by the separate 
extra factors of (/) — 4)a and (Z)— 4)a' combining to give, 

{D-A){a - a) = {D-4:)aa'H At] . (200) 

Of course this makes the sum integrable in D = 4 dimensions, at which point 
we can take D = 4 and the result vanishes on account of the overall factor of 
(D-A). 
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5.4 Sub-Leading Contributions from i6Ac 



The point of this subsection is to compute the contribution from replacing 
the graviton propagator in Table El by its residual C-type part, 



a/3 ^ per 



{D-2){D-3) 



D-3 



D-3 



iSA 



c- 



(201) 

As in the previous sub-sections we first make the requisite contractions and 
then act the derivatives. The result of this first step is summarized in Ta- 
ble El We have sometimes broken the result for a single vertex pair into 
parts because the four different tensors in ()201|1 can make distinct contribu- 
tions, and because distinct contributions also come from breaking up factors 
of 7"J^^. These distinct contributions are labeled by subscripts a, b, c, etc. 

Here i6Ac{x;x') is the residual of the C-type propagator (jU^ after the 
conformal contribution has been subtracted, 



S(f-3)r(^. 




/^2- — 



H^-^ T{D-3) 



(47r)f r(f 
- (n + 1 



T{n+D-3)^yy]^ 



T(n + 



D■^ 



(202) 



As with the contributions from iSAsix; x') considered in the previous sub- 
section, the only way i6Ac{x;x') can give a nonzero contribution in D = A 
dimensions is for it to multiply a singular term. That means only the n = 
term can possibly contribute. Even for the n = term, both derivatives must 
act upon a Ax"^ to make a nonzero contribution in D = 4 dimensions. 

Those of the [2—1] and [2—2] vertex pairs which are not proportional to 
delta functions after the initial contraction of Table El all contrive to give 
delta functions in the end. This happens through the same key identity psip 
which was used to reduce the analogous terms in the previous subsection. In 
each case we have finite constants times different contractions of the following 
tensor function. 



d'Ai6Ac{x;x')dadp 




(203) 
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I 


J 


sub 


tVj'-^ix) i[S]ix; x') tVrix') U^T^^J t5Acix; x') 


2 


1 


a 


(D-3)(D-2) '^'^ ^ ^^^C {x; X) 


2 


1 


b 


-^K^d'{j°do i[S]{x; x')r iSAcix; x')} 


2 


2 


a 


2(D-3)(D-2)'^^^^^(^ x')iSAc(x;x) 


2 


2 


b 


~nn^-\\ ^V^o ^^{x-x') iSAc{x- x) 


2 


2 


c 


+ ^'^'d'^b'"dot[S]ix;x')rtSAcix;x')} 


2 


2 


d 


-UM)^'W9ot[S]{x;x'h't6Ac{x;x')} 


2 


3 


a 


4(r>-3)(D-2) '^^ ^^(^ ^ ^^^c(a;; a;') 


2 


3 


b 


+zr7f^'«^ S^{x-x')Yd'A5Ac{x;x') 


2 


3 


c 




2 


3 


d 




3 


1 


a 




3 


1 


b 




3 


2 


a 


4P-''3k12)^'^;{^^'^^c(x; x') i\S\ {x- a;')7n 


3 


2 


b 


-i(§Ei)'^'5a^^'5Ac(a;; x') t\S\{x- x')^^} 


3 


2 


c 


^^^cix- x') i{S\ {x- x')-f^} 


3 


2 


d 


+i(§^)«29ar9, ibAc{x- x') i{S\{x- x')f} 


3 


3 


a 


siii-3)lL2) [S](x; x')d, ^'i5Ac{x-x') 


3 


3 


b 


-\{%^)K''r^[S]{x■ x')d^iW'^ i6Ac{x; x') 


3 


3 


c 


\{'^-\)K^fi[S]{x- x')d, ^iSAcix; x') 


3 


3 


d 


x')d^^^d'^ t6Ac{x- x') 



Table 16: Contractions from the i5Ac part of the graviton propagator. 
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H^-' T{D-3) ^, 
(47r) 



r(f; 



Ax 



D-2 



H 



D-2 



D 



3 r 



D 



1 1 (aa 



3) 4(D-3) 
1 



, (204) 



IGtt 



Ax2 



16 



Va(^d^S\x-x') + 0{D-A) . 



0(D-4), (205) 
(206) 



It remains to multiply ()206|) by the appropriate prefactors and take the 
appropriate contraction. For example, the [2—1] 5 contribution is, 



K 



D-3 



X 



^ r(f-i) 

, D 

Air 2 



X ^Xl^^' X -^^a/35^5'(x-x') 



TT X -I 

167r2 4 



'(x-x) +0(D-4) . (207) 



We have summarized the results in Table El along with all terms for which 
the initial contractions of Table ITBl produced delta functions. The sum of all 
such terms is, 



€11 



'''^\-y-\^yix-x') + OiD-A) . (208) 



J'"'''^ ' 16772 L 8' 4 
All the lower vertex pairs involve one or more derivatives of i6Ac 

D Ax' 







diiSAc 






3271^ 






doiSAc 






32nT 






d'otSAc 




32nT 



-{D-Q){D-4){aa' 



(D-6)(L>-4)(aa 
{D-6){D-4){aa 



Ax-^-2 



-d[i5Ac , (209) 
aH 



Ax^-2 2Ax^-4 J 

At] a'H 

~Ax^-2~2Ax^-4 



, (210) 
|.(211) 



Their reduction follows the same pattern as in the previous two sub-sections. 
Table ITHl summarizes the results for the case in which all derivatives act upon 
the conformal coordinate separation Ax^. 

When summed, three of the columns of Table ITHl reveal a factor of {D—2) 
which we extract. 



in'^H'^^fD^ fD ^{D-2){D-A){D-6) 



2^7VD 



r 



iD-3) 



aa ) 



r\2- 
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T 
i 


J 


sub 


16^2 90 [X X) 


16,r2 ^ ^ ^ ) 


2 


1 


a 


1 

~2 





2 


-1 

1 


b 


1 
4 


U 


9 


9 


a 


1 
4 


n 

U 


2 


2 


b 


1 

2 


1 
2 


2 


2 


c 


1 

8 





2 


2 


d 


1 
4 


1 
4 


2 


3 


a 








2 


3 


b 








Total 






3 
8 


1 

4 



Table 17: Delta functions from the i6Ac part of the graviton propagator. 



Ax2^- 



We partially integrate (^T^ with the aid of (I17()II171|) and then take D = 4, 
just as we did for the sum of Table HSl 



—I 



x: X 



V(^\vi^ x p-2)(D-4)(D-6) 
V2/ V9 / (D-3) 



/N2-4f 







D2 


7*Axi 


I 4 Ax2^-2 


L VD-2^ 


2(D-1)J 


Ax2^-2 



V2 


/ 7^Ax^ X 


V2 


/ 7*Axi X 


8P-1) 




+ 4(1.-1) 





287r^ ' 



{D-2){D-Q) ^ o\ D{D-l) 



{D-l){D-3y 



8(B-2)2 



D-4 



16(D-2) 



j Ax2^- 



2%4 L2 



,(213) 



(214) 
(215) 
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T 
i 


J 


sub 






||a;|p7'^Ax,j 


||a;|p7*A2;i 










o 
z 


Q 
O 


c 


-{D-lf 


D(D-l) 








2 


3 


d 





(D-l)(D-2) 


D(D-2) 


-2D(D-2) 


3 


1 


a 


4(D-1) 











3 


1 


b 


-2(D-1) 


-2(D-4) 








3 


2 


a 


-2(L)-l) 











3 


2 


b 


2(D-l)(D-2) 


-2(D-l)(D-2) 








3 


2 


c 




(D-4) 








3 


2 


d 





-(2D-3)(D-2) 


-2(D-l)(D-2) 


4(D-l)(D-2) 


3 


3 


a 













3 


3 


b 




ip-l)(D-4) 








3 


3 


c 


i(D-l)2 


i(D-l)(D-4) 








3 


3 


d 


-i(D-l)(D-2) 


(D-2) 


i(D-2)2 


-(D-2)2 


Total 






i(D-l)(D-2) 


2(D-l)-D(D-2) 


-i(D-2)2 


(£)-2)2 



Table 18: iSl^c terms in which all derivatives act upon ^x^{x]x'). All con- 
tributions are multiplied by ^r(f )r(f - 1) (-D-g^^-e) {^aa'f~^ . 
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As already explained, terms for which one or more derivative acts upon a 
scale factor make no contribution in D = 4 dimensions, so this is the final 
nonzero contribution. 



6 Renormalization 

The regulated result we have worked so hard to compute derives from sum- 
ming expressions (HHl), dHHI), (fTT^ - (fTT^ - (HZEl), dUHl) and 
(EIS), 



—I 



+b2H^^ + b3H^^\6^{x-x') + X -?,\n{a)^6\x-x') 



IK 



2%4 



, aa 



Ax2 



ln(yu2Ax^ 



Ax2 



1 



Ax2 



Ax2 

0(D-4). (216) 



The various D-dependent constants in ()216p are, 

.D-A r(f-i) r . 2 ^ 



b2 



287rf 




3)P- 


-4) 




r 


(f + 1) 




2%T 




3)P- 


-4) 




r 


(f + 1) 





-2L) + 1- 



D-2 



D D-1 D-3 



-D + 3 + 



9 



297rf p-3)(D-4) [ ■ - Z}-3 

r(L)-l) f (D + l)(L)-l)(D-4) TT 



(47r)T r(f) 



X — cot 

2{D-3) 2 V 2 / 

{D-l){D^-8D^ + 23D-32) 7' 

48 



,3 . £!l£Ezl)|!,B 



8(Z^-2)2(D-3)2 

TT 



2 X TT /TtDn 



X 



(47r)f r(f) \4V 2 ~-^V 2 

3 (D2-6D + ^ 



4(D-2)2(D-3)2 2j' 



(217) 
(218) 
(219) 



(220) 



(221) 



53 



In obtaining these expressions we have always chosen to convert finite, D = A 
terms with d"^ acting on into delta functions, 

1 



(222) 



All such terms have then been included in 62 and 63. 



Fig. 3: Contribution from counterterms. 

The local divergences in this expression are canceled by the BPHZ coun- 
terterms enumerated at the end of section 3. The generic diagram topology 
is depicted in Fig. 3, and the analytic form is. 



-ictm 



(223) 



1=1 



-K^{ai{aa')-^^d^ + a2D{D-l)H^^ + a3H^^}6^{x-x') .(224) 



In comparing ()216|1 and ()224|1 it would seem that the simplest choice for the 



coefficients is. 



tti = Pi 



0L2 



and as = /^a+^s 



(225) 



D{D-\) 

This choice absorbs all local constants but one is of course left with time 
dependent terms proportional to In(aa'), 

1 In(aa') 



Pi{aa' 

p2{aa'f-^ +b2- D{D-l)a2 



26^2 



0{D-A) , (226) 



^3(aa 



+^\n{aa')+0{D-A) , (227) 
2"7r^ 

1 



"3 



26772 



\n{aa') + 0{D-A) . (228) 



Our final result for the renormalized self-energy is. 



fln(aa' 



15 



^^'^^ 267721 aa 



'^^x — x'] 



IK 



28774 



[aa') 



f\-i 



pln(/i2Aa;2)^ ^ in^H'^ 



Ax2 



28774 



\n{\H^Ax^] 



rln(yu2Ax2 



Ax2 



-7^V' 



Ax^ 



(229) 
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7 Discussion 



We have used dimensional regularization to compute quantum gravitational 
corrections to the fermion self-energy at one loop order in a locally de Sitter 
background. Our regulated result is ()21(jj) . Although Dirac + Einstein is 
not perturbatively renormalizable [3^ we obtained a finite result ()229|1 by 
absorbing the divergences with BPHZ counterterms. 

For this IPI function, and at one loop order, only three counterterms are 
necessary. None of them represents redefinitions of terms in the Lagrangian 
of Dirac + Einstein. Two of the required counterterms (jU^ are generally co- 
ordinate invariant fermion bilinears of dimension six. The third counterterm 
(18 1|) is the only other fermion bilinear of dimension six which respects the 
symmetries fl38ll43|) of our de Sitter noninvariant gauge ^7\i and also obeys 
the refiection property ()80p of the self-energy for massless fermions. 

Although our renormalized result could be changed by altering the finite 
parts of the three BPHZ counterterms, this does not affect its leading behav- 
ior in the far infrared. It is simple to be quantitative about this. Were we to 
make finite shifts Aa^ in our counterterms ()225|) the induced change in the 
renormalized self-energy would be. 



x; X 



aa 



'^^d' + 12Aa2H'^+Aa3H'^\6\x-x') . (230) 



No physical principle seems to fix the Aoj so any result that derives from 
their values is arbitrary. This is why BPHZ renormalization does not yield 
a complete theory. However, at late times (which accesses the far infrared 
because all momenta are redshifted by a(t) = e^*) the local part of the 
renormalized self-energy (j229p is dominated by the large logarithms, 

«2 [InjooO + }^i^^aa')H^^-7\n{aa')H'^]5\x-x') . (231) 



267r2 



The coefficients of these logarithms are finite and completely fixed by our 
calculation. As long as the shifts Aa^ are finite, their impact ()23()j) must 
eventually be dwarfed by the large logarithms ()23H) . 

None of this should seem surprising, although it does with disturbing 
regularity. The comparison we have just made is a standard feature of low 
energy effective field theory and has a very old and distinguished pedigree 
|HlE[linilIIllI3IISllI3II3linilIZllIHllIll2niinil221- Loops of massless 
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particles make finite, nonanalytic contributions which cannot be changed by 
local counterterms and which dominate the far infrared. Further, these effects 
must occur as well, with precisely the same numerical values, in whatever 
fundamental theory ultimately resolves the ultraviolet problem of quantum 
gravity. That is why Feinberg and Sucher got exactly the same long range 
force from the exchange of massless neutrinos using Fermi theory fOl as 
one would get from the Standard Model |TT] . 

So we can use ()229|) reliably in the far infrared. Our motivation for 
undertaking this exercise was to search for a gravitational analogue of what 
Yukawa- coupling a massless, minimally coupled scalar does to massless ferm- 
ions during inflation jT]. Obtaining (j229p completes the first part in that 
program. What remains is to use our result to solve the quantum-corrected 
Dirac equation ((T)). We shall undertake that in a subsequent paper. However, 
it seems clear that the dominant effect must come from the following six 
terms, 



— z 



X ■ 



{x;x') — > ^l^jy In(aa')^- 71n(aa')^|5^(x-x') + ^^g^ 
>^-W)\'^].{sW-m^)\'^^]y (232) 



As adumbrated in the Introduction, these terms are only enhanced by a 
factor of ln(a) relative to the classical part of the Dirac equation (^. That is 
much weaker than the a In (a) enhancement engendered by a Yukawa scalar 
[T] but it can still lead to interesting effects. Note that any such effect will 
be independent of assumptions about the existence and couplings of light 
scalars during inflation. 

We have already commented on the importance of the logarithm terms 
fl232j) . During inflation these infrared logarithms are ubiquitous in quantum 
corrections from massless, minimally coupled scalars and gravitons. A heroic 
recent analysis of scalar-driven inflation at arbitrarily high loop order was 
not able to exclude the possibility that they might even contaminate the 
power spectra of cosmological density perturbations 1^31! The proportional 
correction they make in that case must be small because the logarithms 
would only start to grow at horizon crossing, and must cease growing when 
the mode reenters the horizon after inflation. So the largest enhancement for 
a currently observable mode would be ln(a) < 100. This must be set against 
the tiny loop counting parameter of GH^ < 10~^^. 
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The more significant corrections would be to modes which are still enor- 
mously super-horizon. These are also down by the constant GH^, but the 
time-dependent enhancement factor ln(a) could be arbitrarily big. That is 
what we shall study in our follow-up work. Of course loops of such modes 
can also engender large corrections to effective interactions of low dimen- 
sion. These corrections can grow so large that perturbation theory eventu- 
ally breaks down. Starobinskii has advocated gaining quantitative control 
over this regime by summing the leading infrared logarithms at each order 
|66j . With Yokoyama he has given a complete solution for the case of a 
minimally coupled scalar with arbitrary potential which is a spectator to de 
Sitter inflation [H7j . 

The asymptotic late time effect is small in the simple scalar models for 
which the leading logarithm expansion has been summed. However, it is by 
no means clear that the same must be true for more complicated theories 
that also show infrared logarithms such as quantum gravity jSHl IMl 1^ - 
scalar QED [Ml EH] or Yukawa theory [H Another application of our 
result (j2!29|) is to serve as "data" in checking a leading logarithm formulation 
of Dirac + Einstein during inflation jHEl EHI- Because the fermion does not 
itself engender infrared logarithms it may serve as a spectator for what is 
going on in the gravitational sector. In the leading logarithm limit one could 
obtain an explicit operator expression for the fermion, in terms of its own 
free field and the infrared part of the metric. Although the infrared part of 
the metric would not be known to all orders, it is known to lowest order, and 
that would suffice to compare with one loop results such as ()229|) . This might 
serve as an important intermediate point in the difficult task of generalizing 
Starobinskii's techniques to full blown quantum gravity. 

It is well to close with a comment on accuracy. Although parts of this 
computation are quite intricate we have good confidence that (I229|) is correct 
for three reasons. First, there is the fiat space limit of taking H to zero while 
taking the conformal time to be 77 = —e~^^/H with t held fixed. This checks 
the leading conformal contributions. Our second reason for confidence is the 
fact that all divergences can be absorbed using just the three counterterms 
we have inferred in section 3 on the basis of symmetry. This was by no 
means the case for individual terms; many separate pieces must be added 
to eliminate other divergences. The final check comes from the fact that 
the self-energy of a massless fermion must be odd under interchange of its 
two coordinates. This was again not true for separate contributions, yet it 
emerged when terms were summed. 
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